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We present preliminary results for meson spectral functions at nonzero momentum, 
obtained from quenched lattice QCD simulations at finite temperature using the Maximal 
Entropy Method. Twisted boundary conditions are used to have access to many momenta 
p ~ T. For light quarks, we observe a drastic modification when heating the system from 
below to above T c . In particular, for the vector spectral density we find a nonzero spectral 
weight at all energies. 

Spectral functions provide insight into the real-time dynamics of the quark-gluon plasma. 
In the context of heavy-ion collisions and RHIC, lattice QCD studies carried out so far 
address, among others, the survival of charmonium up to T ~ 1.6T C [ d 121 El E] , the rate 
of dilepton production [E|, and transport coefficients [ 13 IH1 El ED] • Recently, analytical 
studies of spectral functions in strongly coupled gauge theories have also become available, 
using the gauge-gravity duality conjecture for TV = 4 supersymmetric Yang-Mills theory 

[HUE21. 

Meson spectral functions p//(a;,p) can be obtained from euclidean correlators obtained 
in lattice simulations by inverting the relation 

ri \ r duJ is( \ t \ ut \ cosher - 1/2T)] 

Gh{t,P)= —K(r,u)p H (uj,p), K(t,u) = . , 1 

J 2n sinh(u;/2T) 

using the Maximal Entropy Method (MEM) [US]. Here G H (r, x) = (J H (r, x) J ] H (0, 0)) and 
Jh(t, x) = g(r, x)r//g(r, x) with = {1, 75, 7 M , 757^}- Usually only zero momentum 
(p = 0) is considered. In order to better understand the resulting spectral functions, 
especially at high temperature, it may be beneficial to consider spectral functions at 
nonzero momentum as well. There are several reasons for this [E]- At high temperature, 
spectral functions are no longer expected to be sharply peaked. Instead one expects them 
to have a rich structure and contain e.g. a contribution below the lightcone {uj < p = |p|) 
and a momentum-dependent threshold above the lightcone. For very soft momentum 
and energy (u ~ p <^ T) one expects hydrodynamical structure in spectral functions 
related to conserved currents (see refs. [EJEI] f° r more details on both continuum and 
lattice meson spectral functions at infinite temperature; for other recent weak-coupling 
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Figure 1. Euclidean correlator at nonzero momentum normalized with the correlator at 
zero momentum, G(t, p) /G(r, 0) , in the pseudoscalar channel below T c (left) and in the 
scalar and pseudoscalar channel above T c (right). Because of chiral symmetry restoration, 
the latter are hardly distinguishable. The momenta are pL = 2.0, 3.14, 3.72, 4.25, 4.36, 
5.2 (top to bottom). 



calculations, see e.g. ref. [US])- Finding specific momentum-dependent features in spectral 
functions obtained via MEM may help in the interpretation. 

On a lattice with N a sites in a spatial and N T sites in the temporal direction, the 
smallest nonzero momentum in units of the temperature is p/T = 2TTN T /N a . To have 
access to the hydrodynamic regime with many momenta p/T < 1 at high temperature 
(with large enough N T for MEM purposes) is extremely demanding. For that reason 
we have adopted so-called twisted boundary conditions, following closely the formulation 
given in ref. [QZj. Combining the twisting with a standard spatial Fourier transformation, 
the meson momentum reads pL = 27rn — (6\ — 62), where n are integers, 6\ and 62 are 
the twist angles of the two quarks, and L = aN a . In this contribution we present results 
on two lattices, below and above T c m 270 MeV, specified by 

cold: 48 3 x 24, (3 = 6.5, a' 1 ~ 4 GeV, T ~ 160 MeV, 100 propagators, 
hot: 64 3 x 24, (3 = 7.192, a' 1 ~ 10 GeV, T ~ 420 MeV, 100 propagators. 

The estimates for the lattice spacing and the temperature are taken from ref. [12]. Note 
that N T is kept fixed, whereas the other parameters vary. Combining four different twist 
angles in a number of ways, we have ~ 20 different momenta in the range < p/T < 4 
on the hot lattice. The results shown here are for light staggered quarks, am = 0.01 
{m/T = 0.24). 

In fig. [T] the pseudoscalar correlators both in the cold and the hot phase are presented. 
In order to see the relative effect of increasing the momentum we show the ratio of G(r, p) 
and G(t, 0). While in the cold phase the relative effect is large (up to 45% for the largest 
momentum shown), in the hot phase the relative effect is much less (up to 15% for the 
largest momentum shown). In fact, in the hot phase this effect is of the same order as the 
one obtained in a free quark calculation at infinite temperature (note that in the ratio 
the operator renormalization factors cancel). In the hot phase we also present the scalar 
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Figure 2. Spectral functions in the vector channel pu(uj , p) / uj 2 (summed over i 
for various values of p below T c (left) and above T c (right). 



1,2,3) 




Figure 3. Left: vector spectral functions above T c , pu(uj,p)/T 2 , for various values of p. 
The sharp peak near u = is an artefact of the MEM analysis. Right: blow-up of the 
small-energy region with the sharp peak near zero removed. 



correlator to demonstrate chiral symmetry restoration above T c . 

The MEM analysis for staggered fermions is more complicated than for Wilson fermions 
due to the mixing between two signals in the correlators. The equivalent of relation 
for staggered fermions reads 

G s i ag (r,p) = 2 r^Kfau) [ PH (co,p) - (-1)^(^)1 , (2) 
Jo An 

where pu is related to pu by changing Y H to Y H = j^Th. In practice we perform an 
independent MEM analysis on the even and odd time slices, obtaining p e ven — Ph — Ph 
and p dd = Ph + Ph and combine these to get p#. 

In fig. we show spectral functions in the vector channel (Th = 7», summed over 
i = 1, 2, 3), normalized with uj 2 , in the cold and the hot phase. In the cold phase there is 
a clear bound state peak whose position shifts to larger energy with increasing momentum. 
The second broad structure at uj /T ~ 30 (au ~ 1.25) is probably a lattice artefact, due to 
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the finite Brillouin zone [E]- I n the hot phase the bound state appears to have "melted", 
although some structure remains. This requires further study. Note that the horizontal 
axis covers a huge range < u/T < 50. To focus on the region relevant for thermal 
physics, we show a blow-up of pu, now normalized with T 2 , in fig. El A sudden increase is 
visible, which, especially for the larger momenta, appears to occur roughly when u ~ p. 
This is reminiscent of the threshold behaviour observed at weak coupling. The sharp peak 
at uj ~ is an artefact of the MEM analysis and depends e.g. on the resolution along the uj 
axis. Removing this peak by hand and zooming in even further yields the result shown in 
fig. H3 (right). This region is of special interest for hydrodynamics and transport, since the 
electrical conductivity can be defined from the slope of this spectral function at vanishing 
energy via the Kubo formula, o = lim^o Pn(^, 0) / (6a;). We note that the spectral weight 
is nonzero for all uj, which is expected at high temperature. Moreover, for the smallest 
momenta there is a hint of a bump at small uj. Whether this indeed corresponds to the 
expected hydrodynamical structure, and in particular whether the electrical conductivity 
can be (reliably [[7]) extracted from the slope, requires further analysis [ITS]. 
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